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Abstract. After describing in short some problems and methods regarding the smoothness of 
null infinity for isolated systems, I present numerical calculations in which both spatial and null 
infinity can be studied. The reduced conformal field equations based on the conformal Gauss 
gauge allow us in spherical symmetry to calculate numerically the entire Schwarzschild-Kruskal 
spacetime in a smooth way including spacelike, null and timelike infinity and the domain close 
to the singularity. 



1. Motivation 

The main purpose of this article is to study the feasibiUty of the conformal Gauss gauge in 
numerical applications. The importance of this gauge lies in its usefulness in the analysis of 
the gravitational field in a neighbourhood of spatial infinity, as we discuss in the following. For 
review articles see [12l [6] . 

1.1. Isolated Systems 

The motivation to study the regularity of null infinity for asymptotically flat spacetimes arises 
from interest in the following question: How can we describe the gravitational field of isolated 
systems in a rigorous and efficient way? 

Isolated systems in general relativity are self-gravitating models of astrophysical sources. 
Depending on the context, these can be planets, stars, black holes and even galaxy clusters. It 
is not the source that makes a system isolated, but the asymptotic region. The gravitational field 
of isolated systems becomes weak as we move away from the source. As nothing is completely 
isolated from the rest of the universe and as we can not yet test our assumptions by experiment, 
details of this idealization need to be decided upon in accordance with our physical intuition, 
expectation and most importantly, in accordance with the equations describing the model. 
Isolated systems are modeled in general relativity by asymptotically flat spacetimes in which 
the metric on the complement of a spatially compact set approaches near infinity the flat metric 
in a specified way. 

Many physical concepts like mass, linear and angular momentum are defined in the asymptotic 
region with respect to infinity. The use of coordinate dependent limits to infinity for calculating 
these quantities can deliver misleading results if one is not careful enough. A geometric, 
coordinate independent way to deal with isolated systems has been proposed by Penrose [15] by 
use of conformally rescaled spacetimes with a smooth conformal boundary. The main idea is to 
adjoin infinity to the physical spacetime by a conformal rescaling of the metric, as is done in 
complex analysis, where the Riemann sphere is constructed by adjoining infinity to the complex 




Figure 1. Numerically generated Penrose diagram of the Minkowski spacetime with Cauchy 
and hyperboloidal surfaces 



plane [16]. One studies the equivalence class of the physical metric g given by the conformal 
rescaling g = O'^g with a function ^ > 0. As the physical distances go to infinity the conformal 
factor approaches zero at a specific rate so that the rescaled metric g is regular at the set 
{9 = 0} of the conformal extension if appropriate fall-off conditions are satisfied. This geometric 
construction replaces coordinate dependent limits to infinity by local differential geometry at 
{0 = 0} and allows an elegant description of the asymptotic structure of isolated, self- gravitating 
systems which has proven very fruitful in mathematical analysis. 

The Penrose diagram in Fig. [1] depicts the causal structure of the Minkowski spacetime 
illustrating the different asymptotic regions. One can go to infinity along timelike, null or 
spacelike directions approaching timelike i^, null or spatial infinity . In such a Penrose 
diagram, light rays are straight lines with 45 degrees to the horizontal. Plotted are also Cauchy 
surfaces that approach spatial infinity and hyperboloidal surfaces that cut null infinity. The 
conformal compactification of the Minkowski spacetime, as found for example in [6j, is the 
simplest application of the idea of conformal compactification. 

1.2. Spatial infinity 

The details of the conformal compactification technique developed by Penrose have been 
motivated rather by explicit examples or well studied assumptions about the fall-off behaviour 
of certain fields than a detailed study of the full non-linear Einstein's field equations, which 
poses the following question [T7]: How general is the description proposed by Penrose? 

We would like to know by some general argument whether we have a large class of non-trivial, 
asymptotically fiat, radiative spacetimes that admit a smooth conformal boundary such that we 
can apply the conformal compactification technique. To answer this question, the solution space 
to the Einstein equations needs to be studied with an emphasis on the asymptotic structure of 
the gravitational fields. 

The available explicit solutions are not general enough to study the solution space by direct 
means, therefore we are led to abstract analysis. By sufficient knowledge of the properties of 
the spacetime, we can get some general results. For example we know that asymptotically fiat, 
vacuum, stationary spacetimes admit an analytic compactification at null infinity [5j. For more 
general results that include radiative spacetimes, the initial value problem needs to be studied. 

Unfortunately, the analysis of the initial value problem for the Einstein equations in the 



compactified picture is not straightforward, as the Einstein's equations are not conformany 
invariant and compactification leads to formany singular equations. However, the equations are 
conformally regular as Friedrich showed by constructing a system which is equivalent to the 
Einstein's equations for 6 > and is regular for all values of the conformal factor so that the 
equations can be analysed on the conformal extension of the spacetime [7j. With this system 
Friedrich showed that if conformal data on a spacelike hypersurface that cuts null infinity is given 
such that it is smooth up to null infinity, then the evolution of this data has a smooth conformal 
boundary [8j. This analysis showed that the decision on the smoothness of null infinity is made 
at spatial infinity. 

By the gluing techniques developed by Corvino and Schoen [2] one knows that there exists 
a large class of non-trivial initial data which is Schwarzschild or Kerr in a neighbourhood of 
spatial infinity. Combined with Friedrich's earlier results this leads to the existence of a large 
class of non-trivial radiative vacuum spacetimes with smooth null infinity [3] . 

However, these spacetimes have a special asymptotic structure. The question still remains 
whether more general spacetimes exist. One would also like to construct such spacetimes 
numerically. These tasks require a more sophisticated method to study spatial infinity. 

1.3. Reduced Conformal Field Equations 

A detailed study of the solutions in a neighbourhood of spatial infinity is complicated by the 
fact that in general point compactification at spatial infinity leads to singular conformal data. 
However, Friedrich was able to pose a regular finite initial value problem near spatial infinity 
by using the reduced conformal field equations that he introduced in [9]. These equations are 
based on the conformal Gauss gauge in which spatial infinity in an asymptotically flat spacetime 
can be represented as a cylinder [10]. Subsequent analysis showed that in general, logarithmic 
singularities arise in a small neighbourhood of spatial infinity. Friedrich obtained necessary 
regularity conditions on the Bach tensor. Later on, Valiente Kroon showed that these conditions 
are not sufficient and obtained further obstructions to smoothness of null infinity [IBJ • 

While the question about the necessary and sufficient conditions for smoothness of null infinity 
still remains, we can ask whether one can deal with some mild singular behaviour numerically by 
using the reduced conformal field equations, which turned out to be the proper tool to study the 
asymptotic behaviour of solutions in analytic work. A strong interaction between mathematical 
and numerical studies in this question might give new impulses in both directions. The desire to 
have numerical simulations which allow one to understand and control the asymptotic structure 
of spacetimes is the main motivation behind my work. 

In the following, we will summarize ingredients of the underlying conformal Gauss gauge for 
the reduced conformal field equations. 

2. Conformal Gauss gauge 

The main reference for the properties of the conformal Gauss gauge that we discuss below is 
This gauge is based on conformal geodesies, which are autoparallel curves with respect to 
a Weyl connection. 

2.1. The Weyl connection 

The main step from Euclidean geometry to Riemannian geometry is the removal of the 
assumption of integrability of vectors by parallel transport. A vector parallel transported along 
a closed curve changes in general its direction due to curvature, while its norm stays the same. 
Weyl, arguing that the choice of the unit of measurement, i.e. the gauge, should also be subject 
to local variations, demanded "the non-integrability of the transference of distances" [19j. This 
simple requirement leads to a rich geometry, which Weyl used in his attempt to unify the 



gravitational and electromagnetic fields. While his attempt was not successful, the main idea 
has been the basis during the construction of modern gauge theories |14j . 

Driven by such conceptual and also other philosophical considerations, Weyl introduced a 
torsion free connection which is not necessarily the Levi-Civita connection of a metric but 
preserves the conformal structure and thus satisfies 

^aQfiu = -"^faQfiu- (1) 

If the 1-form / is exact, it can be written as / = Q^^dQ, and the Weyl connection V will be the 
Levi-Civita connection of a metric in the conformal class 5^,^ = Vt^g^y. The Weyl connection 
V is invariant under a conformal rescaling g = 'd'^g with a function > 0, so the 1-form / 
transforms according to / = / — di?. 

The relation of the Weyl connection V with a metric connection V of a metric g in the 
conformal class is given by 

V - V = where ^ 5" J, + 5",!^ - g^^g'^fx. 

We see that the Weyl connections are characterized by 1-forms /. Contrary to Weyl who used 
the additional freedom introduced by the 1-form / to construct a unified field theory, Friedrich 
exploited it as a gauge freedom of the conformal structure in his construction of the reduced 
conformal field equations, which are equivalent to the Einstein's field equations. 

A torsionfree connection implies the notion of parallel transport and allows the construction 
of geodesies, which we will study next. 

2.2. Conformal geodesies 

Null geodesic congruences, when they are smooth, provide a valuable tool to study the 
asymptotic and causal structure of spacetimes. As they are invariants of the conformal structure, 
one might assume that time- or spacelike curves that are conformal invariants might also be 
useful in such studies, which is indeed the case. Conformal geodesies are conformally invariant 
in the sense that they are independent of the metric chosen in the conformal class. 

A solution to the conformal geodesic equations does not only provide a spacetime curve, but 
along the curve also a Weyl connection, a conformal factor and a frame which is orthonormal 
for a metric in the conformal class. While the equations are independent of coordinates, we will 
have to write them in some coordinate system, as we will study them in numerical applications. 

Given a metric g, the equations that define the conformal geodesic x^{t) are written for its 
tangent vector x^[t) and a covector /^(r) that defines a Weyl connection 

(Vixr + S(/)/^xV = 0, 

where L^j^ = ^ Rij — ^ Rgij is the Schouten tensor. Written in terms of the Weyl connection 
that is defined by the 1-form, V = V + S{f), the first equation reads V^i = 0. The equation 
for the spacelike frame fields eJ^{T) along the curve is V^Ca = 0. 

From ([1]) we see that on a given curve x^(r), the metric O'^g^^u with a conformal factor > is 
parallel transported with respect to V if satisfies V ±0 = f^x^^ along the curve. A remarkable 
property of the conformal geodesies is that in a vacuum spacetime the equation for the conformal 
factor 6{t) can be solved explicitly, such that 0{t) is known a priori in terms of initial data along 
the conformal geodesies. For certain initial data, the conformal factor can be written as 



is a conformal factor on the initial hypersurface, hij is the with 17 rescaled spatial metric in 
compactifying coordinates, k is a free function determining the time coordinate at which the 
conformal geodesic cuts J^. One can also show that 

d,:=9Ue,^=(-2T^,-eam). (3) 

For a derivation of these results, see [S]. 

Some of the conformal geodesies leave the spacetime through null infinity. This might seem 
to be a counter-intuitive behaviour as the curves are timelike everywhere. The explanation is 
that they experience a non- vanishing acceleration by the 1-form /. A similar behaviour show 
hyperboloidal slices which are spacelike everywhere but also reach null infinity. One should 
remember that the causal nature of a curve or a surface does not completely determine its 
asymptotic behaviour. 

To construct a conformal Gauss gauge one uses conformal geodesies in a similar way as one 
uses metric geodesies to construct the Gauss gauge. One specifies a congruence of timelike 
vectors, a conformal factor and a 1-form on an initial hypersurface. The timelike geodesies 
starting from this surface provide the conformal Gauss coordinates. Spatial coordinates are 
dragged along. In [11] Friedrich constructs conformal Gauss coordinates in the Schwarzschild- 
Kruskal spacetime and shows that they do not only cover the global solution, but also the 
conformal extension in a smooth way. We will construct this gauge numerically. 

2.3. Numerical experiments 

First we describe the general route to construct the gauge on a given background and then we 
mention the numerical results. 

Assume a solution to the Einstein field equations (A^,^) has been given. A procedure to 
construct a conformal Gauss gauge with timelike conformal geodesies and an orthonormal frame 
along them ( x^(r), /^(r), e^(r) ) for a given spacetime solution to the vacuum Einstein equations 
can be given as follows: 

(i) Find a conformal compactification of the global solution g = 'd'^g. In general, you will need 
different coordinates in different asymptotic regions. 

(ii) On a spatial slice with a spatial metric /i, introduce compactifying coordinates and rescale 
h with an appropriate function fi, such that in these coordinates you have h = Q^/i. The 
metric h is needed for the calculation of the conformal factor ([2]). 

(iii) Set initial data (x''(0), ±'^(0), /^(O), e^(0)) on the initial hypersurface r = according to 

/^(O) = Q-^d^n, ^(O)x'^(O) = 0, 9''~g{x,x)\r=o = %~g{x,x)\r=o = -1, 

6'2g(ea,x)|r=o = 0, 6'^5(ea,ea)|r=o = 1, with a = 1,2, 3. 

The timelike frame vector is given by x itself. The frame is not unique, the freedom 
corresponds to the freedom of Lorentz rotations. One also needs to choose the free function 
K that determines the form of in the conformal Gauss gauge. 

(iv) Solve the following system of ordinary differential equations with the above initial data 

{drxY = 

[d^xY = -r/^i^i^ - 2{f,x^) x'' + {gxpi^xn 5^^/,, 

{drCar = -r/.i'e,^ - (/.e^) - iMne>: + {gxpe^i^ g'^'^f,. (4) 



(v) Check the quahty of the solution using ^ and ([3]) . 

The right hand sides of ([H) are calculated using the computer algebra package MathTensor. For 
the integration of the ODE system, I used a 3th order Runge-Kutta integration algorithm. 

2.3.1. Schwarzschild-Kruskal We calculate the conformal Gauss initial data for the 
Schwarzschild-Kruskal spacetime. The resulting gauge is illustrated in a numerically generated 
conformal diagram Fig. [2j The physical Schwarzschild metric for r > 2m reads 

5 = - (l - ^) dt^ + (l - ^) ^ dr' + da'. 

We transform the metric using the retarded and advanced null coordinates 

w = t — {r + 2m In (r — 2m)), v = t + [r + 2m In (r — 2m)). (5) 

After the inversion z = ^ and the rescaling with "d = z following g = ■d'^g we get 

g = —2^(1 — 2mz) dw'^ + 2 dw dz + da', g = —z'(l — 2mz) dv' — 2 dv dz + da' 

The coordinates w,v extend analytically into regions where r < 2m. For a simulation through 
we use the retarded null coordinate w, for a simulation through the future horizon to the 
future singularity we use the advanced null coordinate v. 

We give initial data on a t =const. slice. From the transformation ([5]) we have 

2^(1 — 2mz) ' z'^{l — 2mz) ' 

which leads to the induced spatial metric 

^=.4(1 .'2^,) dz' + ^da'. 

From now on, we choose m = 2 and present the calculation in the outgoing case to simplify 
the formulae. We compactify h using the conformal factor Vl = , so that h = ^'h \s 



1-22' 

diffeomorph to the standard metric on the three sphere h = dx' + sin^ xda', as can be shown 

sin X 
2(l+sinx)- 



by the transformation z{x) ~ — 



The initial data for the 1-form / can be given as /(O) = z(i-2z) d^' orhogonality condition 
fij.x^ = and the normalisation requirement delivers i;(0) = ^^^^==k(9^. The timelike frame 
vector is given by eo = x. Initial data for the spatial frame vector reads 

1-2^. ^ (l-2zVr^ 
= ' + ' 2 

We can also take these steps using other compactifications of Schwarzschild, for example the 
isotropic compactification. The use of different coordinate systems or different compactifications 
for calculating the right hand side of ^ does not effect the quality of the calculation. The 
resulting conformal diagram depends on the choice of the spatial coordinate and the free function 
K. Fig. [2] is taken from [13] and is the result of a numerical calculation where k has been chosen 
such that J^"*" is a straight line in the corresponding conformal Gauss gauge. Illustrated is 
the "upper right part" of the Penrose diagram for the Schwarzschild-Kruskal spacetime. The 
lower horizontal line corresponds to the hypersurface \t = 0} in the standard Schwarzschild 




coordinates, where also {r = 0}. We see that the conformal geodesies cover in a smooth way 
spacehke, null and timelike infinity and the domain close to the singularity. 

A typical question regarding calculations going up to infinity is whether it is not "too far 
away". After all, the Earth where we are doing our measurements is not infinitely far away 
from the sources and we do not move along null geodesies as an observer along null infinity 
does. In the Schwarzschild case, we have a length scale M at our disposal, so we can compare 
the astrophysical situation with the numerical calculations. In the conformal diagram, three 
curves with constant distances to the event horizon d = 2M, lOM, 300M have been plotted. 
Numerical codes today have outer boundaries at about 300M to lOOOM. The distance lOOOM 
corresponds to about 3000km for a two solar mass black hole. Even if we consider supermassive 
black holes, the corresponding distance is incomparably small with respect to the thousands 
or millions of light years that separates us from the astrophysical sources. As one sees in the 
Fig. El already the curve at 300M is barely distinguishable from J^+. In this representation, 
the numerical effort for simulating the region from 300M to J^^ is very small whereas in the 
standard approach putting the outer boundary from 300M to lOOOM costs considerable effort 
in terms of numerical techniques and computational sources. Therefore, one can conclude that 
the conformal compactification technique is more apt to deal with the asymptotic region. 



2.3.2. Kerr Going beyond spherical symmetry, we can solve the conformal geodesic equations 
in Kerr spacetime. We take the Kerr metric in Boyer-Lindquist coordinates, do the inversion 
and make an Eddington-Finkelstein-like transformation using ingoing and outgoing light rays as 
coordinate lines as before. For the outgoing case, the metric becomes 

2m\ 2 2 4am . 2/, , 
1 — I aw ^-atf az — sm aw a<p 



+ '^sin^edzd(t) + ^de'^ + ^[[-^ + a^] - Aa^ sin^^ ) sin^^ #^ 




2 



where 



^ 1 2 2 /I A 1 2m 2 
1^ = + a cos 0, A = h a . 

Z^ 2^ z 



Consider the coordinate transformation ziV) = ^ , '^'"'^ — with 6 = \Jw? — c?. Following 141, a 
t =const. hypersurface in the original Boyer-Lindquist coordinates can be compactified with the 
conformal factor 

sinx 8z^ 



(1 — mz)\/l + a'^z'^cos'^ 



U ^ 2 , ■ 2 ,/)2 , ■ 2 ■ 2q(. , 2l+2W(Sf(x)) -2 -2 

h = ax + sm + sm x sm p 1 + a sm x sm 



The compactified spatial metric reads 

S sin X 

We get for the initial data from / = 17^^ dVL 

2 — mz + a? cos^ 6 ^ ^ a? z"^ cos 9 sinO 
z{l — mz) {1 + a^z-^ cos^ 9) l + a^z^cos^^ ' 

and for x 

, , (\-rrLz){\^c?z^co'^9) 
i(0) = \ , \ ^ =Kd^. 
bz^\/\ + arz^ cos^9 — 2mz 

Note that for a = we get the formulae for the Schwarzschild case. 

The conformal geodesic equations @ can be solved numerically with the above initial data. 

The behaviour of the curves is similar to the Schwarzschild case, with the notable difference 

that Cauchy horizons exist. Conformal geodesies pass through the horizons where one needs to 

change the coordinate representation of the metric to integrate the equations further. It would 

be interesting to see effects of this behaviour in a Cauchy evolution. 



3. Reduced conformal field equations 

The equations are written for the following variables: a frame field , a Weyl connection Tj^^-, 
the Schouten tensor Lij = \R{ij) — ~ 9ij^ ^^'^ ^he rescaled conformal Weyl tensor 

W)ki='eC',M- They are 

9t^i^j = -^i^j^io + g^o^ij + g^jLoi - gjoLj^'^ + 9W''jQi, 
dtLij = -T^^Lkj + diW''jQi, 
ViW%, = 0. 

The functions 9 and d^ are known a priori in terms of initial data. Main advantages of this 
system, especially for numerical work, are the following: 

• The system consists mainly of ordinary differential equations except the Bianchi equation, 
which has symmetric hyperbolic reductions. This is advantageous especially when one is 
dealing with complicated geometries. 

• The location of the conformal boundary is known explicitly. 

One difficulty is that the equations become degenerate at the set X"*" where null infinity meets 
spatial infinity [10]. A first step in numerical work with these equations is the spherically 
symmetric case. We know that the conformal Gauss gauge covers the complete Schwarzschild- 
Kruskal solution and is robust enough to do numerical calculations. Therefore, the Cauchy 
problem in spherical symmetry with Schwarzschild initial data can be expected to be solvable 
without difficulties. 



3.1. Spherical symmetry 

To determine the only non-vanishing Weyl tensor component in spherical symmetry, we need to 
solve the following subsystem of the reduced conformal field equations 

drT,\ = -(T,\f + L22~\9W%,, 

drW% = -3T2\W%,. (6) 



The quantities 9 and do are given by 



=- , dor=-2T— . 

The degeneracy of the equations at X+ is not present in the simple case of spherical symmetry. 
The reduced conformal field equations become a system of ordinary differential equations. 

The initial data can be calculated as in [10]. We choose a frame such that h{ea,eb) = Sab 
with respect to h = Q'^h and then rescale it via = KeJ^. The connection coefficients are 
calculated with respect to ej^. For the Schouten tensor and the electric part of the Weyl tensor 
we calculate 

Lij = tij + —f hij^ , W%j = (^^ tij + - Si 

with tij = DiDjQ, — I Sij DkDiQ 5^^ and Sij = fjj — \f5ij are calculated with respect to ej^. By 
solving the system ([6]) we can generate the entire Schwarzschild-Kruskal solution including the 
conformal extension. 

To check the quality of the solution, we compare the Weyl tensor in the frame representation 
adapted to the conformal geodesies. In the conformal Gauss gauge, we calculate the timelike 
coframe using the orthonormality conditions 

pi p 



^0 plpO_plpO' ^1 p ip _ p Ip 

and compare the rescaled Weyl tensor in the conformal Gauss gauge on the Schwarzschild- 
Kruskal background given by 



W%, = a\e,^e,W \c%^p, (7) 



with the result of the numerical simulation of the Schwarzschild initial data. The calculations 
agree with each other within the numerical order of integration, so the resulting spacetime can 
be illustrated by the same conformal diagram Fig. [2] 

While numerical tests of tetrad formulations in spherical symmetry can be quite delicate and 
instructive p!] , the simplicity of ([6j) does not allow us to draw representative conclusions for the 
general case. Still, the study suggests that global numerical calculations for spacetimes with less 
symmetry might be possible using the reduced conformal field equations. 



4. Conclusions 

Basic motivations for numerical studies of conformal compactification are to avoid artificial 
outer boundaries and to have rigorous analysis tools for numerically generated spacetimes with 
efficient codes. A detailed understanding and numerical control of the conformal structure at 
spatial infinity is required for a successful implementation of conformal compactification. The 
tool for this task, i.e. the reduced conformal field equations, has been developed by Friedrich. 
Within this article we have studied the numerical feasibility of the underlying conformal Gauss 
gauge in simple cases. 

The results of this study can be summarized as follows: One can numerically reproduce 
Friedrich's covering of the complete Schwarzschild-Kruskal solution using the conformal Gauss 
gauge. Numerical simulations suggest that one can also cover the complete Kerr solution. For 
the first time, the Cauchy problem for the simplest asymptotically fiat black hole spacetime has 
been solved globally including spacelike, null and timelike infinity and the domain close to the 
singularity. This was made possible by the reduced conformal field equations. 



As the studied examples are very special, they are not representative for the general case. 
One should study the robustness of the conformal Gauss gauge in radiative spacetimes. My 
current work is directed towards numerical simulation of solutions to the reduced conformal 
field equations with a non-vanishing radiation field along null infinity. The imposed geometry 
by the cylinder at spatial infinity and the degeneracy of the equations at requires highly 
developed numerical techniques for the solution of this problem. 
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